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We propose a method for the determination of the interaction potential of Rydberg atoms. Specif- 
ically, we consider a laser-driven Rydberg gas confined in a one-dimensional lattice and demonstrate 
that the Rydberg atom number after a laser excitation cycle as a function of the laser detuning pro- 
vides a measure for the Rydberg interaction coefficient. With the lattice spacing precisely known, 
the proposed scheme only relies on the measurement of the number of Rydberg atoms and thus 
circumvents the necessity to map the interaction potential by varying the interparticle separation. 



Among many fascinating systems encountered in mod- 
ern ultracold atomic and molecular physics are Ryd- 
berg atoms, i.e., highly excited atoms with large prin- 
cipal quantum number n, whose size can easily exceed 
that of ground state atoms by several orders of mag- 
nitude fD. The enormous displacement of the atomic 
charges makes Rydberg atoms highly susceptible to ex- 
ternal fields and at the same time is the origin of their 
strong mutual interaction. In ultracold gases, the strong 
dipole-dipole interaction has been shown theoretically [2] 
and experimentally [T to entail a blockade mechanism in 
the laser excitation of Rydberg atoms, thereby effectuat- 
ing a collective excitation process of Rydberg atoms [4] . 
Two recent experiments even demonstrated the block- 
ade between two single atoms a few micrometers apart 
[5] . The strong dipole-dipole interaction renders Rydberg 
atoms also promising candidates for the implementation 
of protocols realizing quantum gates or efficient multi- 
particle entanglement [6]. In fact, only very recently 
a CNOT gate between two individually addressed neu- 
tral atoms and the generation of entanglement has been 
demonstrated experimentally by employing the Rydberg 
blockade mechanism [7 . Collective excitation of Rydberg 
atoms is currently being studied extensively [8, 9 , pre- 
dicting amongst others the dynamical creation of crys- 
talline photonic states [10 , the formation of fermionic 
collective excitations in a lattice gas of Rydberg atoms 
[TT] , or the ability to act as a universal quantum simula- 
tor [E]. 

The key ingredient for many of the predicted effects in 
ultracold Rydberg physics is the strong and long-ranged 
mutual interaction of Rydberg atoms which has predom- 
inantly van der Waals character at large separations [6]. 
Hence, the precise knowledge of the corresponding in- 
teraction potentials is crucial for making reliable predic- 
tions as well as for correctly interpreting experimental 
data. Accordingly, the long-range Rydberg interaction 
has been discussed in several recent works and theoreti- 
cal predictions were made [HI |T4 . Experimentally, the 
mechanical effect of van der Waals interactions has been 




FIG. 1. Schematics of the considered setup. Ultracold 
atoms in a lattice are excited from the ground state \g) 
(blue/dark gray) to their Rydberg state |e) (red/light gray). 
The Rydberg-Rydberg interaction shifts the excitation from 
its single atom resonance (top of the figure). This dipole 
blockade effect can be compensated by the laser detuning A. 
In this particular example the excitation is resonant to the 
many-body state \geegeeg). 



observed in an ultracold Rydberg gas [15]. 

In this letter we propose a scheme to determine the 
interaction strength between two Rydberg atoms by 
means of the Rydberg blockade mechanism. Specifically, 
we consider the coherent Rydberg excitation of a one- 
dimensional lattice of trapped ground state atoms. A 
thorough analysis of the excitation spectrum allows us 
to identify distinctive values for the detuning of the exci- 
tation laser from the single atom resonance that entail a 
well defined number of Rydberg atoms within the lattice. 
From these detunings the interaction strength of two Ry- 
dberg atoms can be deduced. The possibility to deter- 
mine the interaction strength of two Rydberg atoms by 
investigating laser-driven samples of ground state atoms 
was first mentioned in Ref. [16 where the anti-blockade 
of Rydberg excitation is investigated. Experimentally, 
this anti-blockade effect was observed in [17]. However, 
the experiment was conducted in a frozen gas rather than 
in a lattice and was thus not aimed to predict the inter- 
action strength of Rydberg atoms. 

As in Ref. [18] , we consider an ultracold gas of ground 
state atoms in a one-dimensional lattice consisting of N 
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sites with regular spacing a and a site-independent fill- 
ing of A^o atoms. The lattice is subject to a laser exci- 
tation to the Rydberg state with a single atom Rabi fre- 
quency Qq. While this laser transition is usually of two- 
photon character, the involved intermediate state can be 
removed from the state space by an adiabatic elimination 
if the excitation lasers are strongly detuned with respect 
to this level [19 . Moreover, we assume that the dipole 
blockade is effective within each lattice site k such that 
our system can be described by means of an ensemble of 
individual two- level systems, each of which possessing a 
ground state |^)^ and an excited state | e)^ that are laser- 
coupled with a collective Rabi frequency ft = v^TVo^^o- 
Often, these particular two-level systems where the Ry- 
dberg excitation is shared among Nq atoms are referred 
to as superatoms. A sketch of our setup is provided in 
Fig. [l] Applying the rotating wave approximation, one 
obtains the Hamiltonian 18 



N 



N 



a. 



ik) 



V 



N (k) (0 

l>k ' 



k\' 



(1) 



The operators (J^^\i G {x, ^, z} act only on the super- 
atom located at site k and take on the usual Pauli- 
matrix form when expressed in the local superatom ba- 
sis The excitation number operator reads 

then n'e^ = \ [(J^z^^ + 1]. Hamiltonian ([l]) is valid within 
the frozen gas limit where no movement of the Rydberg 
atoms is assumed. For n = 70 and an initial separa- 
tion of a = 10 /im, the change of the distance of two 
Rydberg atoms due to their repulsive interaction is esti- 
mated to be less than a fraction of 10~^ for an evolution 
time of 1 /is [flO] ; this displacement increases for smaller 
initial distances and higher principal quantum numbers 
(we consider ^''Rb atoms as an ubiquitous example when- 
ever actual numbers are provided). To actively prevent 
vast movement of the Rydberg atoms, one might leave 
the optical trap on during Rydberg excitation. For nSi/2 
Rydberg states with n > 40, ionization lifetimes of Ry- 
dberg atoms due to far infrared trapping lasers are esti- 
mated to reach the /is regime [20 , which is sufficient for 
the present purpose. 

Hamiltonian ([T]) consists of three different contribu- 
tions. The laser coupling of each single atom's ground 
state to the excited state is given by ^^Xl^i^a;^^- 
^A^^-^ai^^ provides the energy gap A between the 
excited state and the ground state of the superatom 
at each individual site. The interaction strength be- 
tween two Rydberg atoms at neighboring sites is given 
hy V = Cm/a^ and is considered to be a perturbation of 
the atomic energy levels. Here, we focus on the regime 
of a dominant van der Waals interaction (m = 6); for 
lattice spacings in the /im regime, the next order contri- 
bution (m = 8) is more than two orders of magnitude 
smaller. Moreover, we assume a repulsive interaction. 



y > 0, which is common for Rydberg atoms in their 
nS 1/2 state JTS . 

For a weak laser coupling, the contributions diagonal 
in the spin basis dominate the Hamiltonian. In this case 
the laser coupling leads to a small off-diagonal pertur- 
bation. We thus assume \Q\ <C V^\Q\ <C |A| and divide 
the Hamiltonian into two parts, H = Hq + H\ group- 
ing together the laser detuning with the next-neighbor 
Rydberg interactions to give the dominant contribution 
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while the perturbation consists of the laser coupling as 
well as the long-range Rydberg interactions 
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Introducing operators for the total excitation number, 
Ne = Yl^=i^^e \ and the next-neighbor excitation pair 
number, N^e = Yl^=i ni^^ni^~^^\ we can rewrite the un- 
perturbed Hamiltonian as 



Ho = A{N,-N/2) + VN,,. 



(4) 



For a given canonical product state \a) = \ s1S2 . . . sat), 
^ {gk^^k} the unperturbed energy eigenvalue is given 

by 



E{a) = A[iVe(a) - N/2] + VN^^ia), 



(5) 



where Ne{a) and Nee{a) denote the integral eigenval- 
ues of the operators and for the state re- 
spectively. In general, there are multiple canonical prod- 
uct states of equal {Ne^Nee) which are degenerate with 
respect to the unperturbed Hamiltonian. Furthermore, 
for specific values of the ratio of A and the simple 
form of Eq. ^ allows for additional degeneracies be- 
tween states belonging to different (A/'g, A'ee)- Figure [2|c) 
presents the energy eigenvalues of the full Hamiltonian 
as a function of a varying ratio A/V. The points of 
high degeneracy only occur when A/V takes on cer- 
tain rational values. For A = the laser is tuned to 
atomic resonance. In this case, for i^o, the canonical 
ground state \ G) = \gg . . . g) is degenerate with all states 
that lack neighboring excitations, i.e., dipole-blockaded 
states. For non-zero detunings, the canonical ground 
state can also be brought to degeneracy with other states. 
Physically, this corresponds to a situation in which the 
lasers resonantly couple \G) to states containing neigh- 
boring excitations, compensating the Rydberg- Rydberg 
interaction by means of the detuning. An example is 
given by A/V = —1/2: In this case \G) is degenerate 
with states containing exactly two neighboring excita- 
tions {\eeg...g) ,\geeg...g) ,...\g...gee)}. The high- 
est degeneracy, i.e., the minimal number of energy bands, 
is achieved for A/V = —1. 
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FIG. 2. (a) Total number of Rydberg atoms and (b) number 
of neighboring Rydberg excitations after the laser excitation 
of a = 8 lattice of ground state atoms as a function of 
the detuning from the single atom resonance; a laser coupling 
strength of Q = 0.151/ is considered. The thin gray lines 
correspond to excitation durations of T = 15, 18, 21, 24, and 
27n~\ The thick red line is an average over many different 
laser cycle times, (c) Energy spectrum for the same setup. 
The dashed line represents the canonical ground state | G) for 
r2 = 0. The peaks of panels (a) and (b) occur at the points 
where | G) crosses many-body states possessing different num- 
ber of Rydberg excitations, as described by Eq. (|6]). 



Here, we are interested in the specific ratios A/F where 
the state | G) becomes degenerate with states possessing 
different numbers of Rydberg excitations. Particularly, 
this occurs at 
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cf. Eq. ([5|. At these detunings, the canonical ground 
state is resonant with the states that satisfy N^e/Ne = 
{k — 1)//^, e.g., all states containing = tz e N>2 Ryd- 
berg excitations whereof N^e = Hi — 1 reside at neighbor- 
ing lattice sites. Performing the Rydberg excitation at 
such defined detunings will thus drive the lattice in states 
that contain an increasing number of Rydberg excitation 
and -pairs. Starting from the canonical ground state at 
t = 0, we solved the time-dependent Schrodinger equa- 
tion associated with Hamiltonian Q. Figure |2|a) shows 
the Rydberg excitation number Ne after the excitation 
cycle {t = T) as a function of the laser detuning from the 



single atom resonance; different durations T of the laser 
excitation are considered. As expected from Eq. (|6|, we 
find at the well-defined values A^^ sharp peaks in the Ry- 
dberg atom number N^. The same holds for the number 
Nee of Rydberg pairs as illustrated in Fig. |2|b). From 
Fig. [2|c) it becomes evident that the canonical ground 
state crosses submanifolds of many-body states at values 
for the detuning other than specified by Eq. (|6|. How- 
ever, these states are only weakly coupled to |G), giving 
rise to very narrow avoided crossings. As a result, unreal- 
istically long propagation times are necessary to actually 
populate these many-body states. Correspondingly, no 
resonance peaks are visible in Figs. |2|a) and (b) at the 
corresponding detunings. 

Relating the detunings A^^ to our theoretical predic- 
tions allows the determination of the van der Waals in- 
teraction coefficient via the relation 
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Since the lattice spacing a is usually known very precisely, 
any uncertainty in the determination of Cq via Eq. ^ 
stems solely from the measurement of the resonant de- 
tunings A^. Because of the rather broad A = peak 
that serves as reference for the detunings, an accurate 
determination of the A/V < resonance positions can 
be difficult. As a remedy, employing the relative posi- 
tion A^^/^+i = A^ — A/^+i of neighboring peaks provides 
a more reliable source for precision measurements. One 
finds Cq = ^{n + l)A^^^+ia^. If unknown, the actual 
value of n can be determined by relating two neighbor- 
ing peaks, A/^/A^+i = 1 — k''^. Experimentally, Ne can 
be determined by field ionization of the Rydberg atoms 
and their subsequent detection via micro-channel plates. 
For the measurement of Nee^ a spatially resolved detec- 
tion is necessary. 

The magnitude of the A^^ - and therefore the exper- 
imental feasibility - is determined by the lattice con- 
stant a as well as the principal Rydberg quantum num- 
ber n. The expected relevant detunings A^^ can be es- 
timated by employing the theoretically determined Cq 
coefficient of Ref. [13 . Taking a = 10 /im as a typi- 
cal example for micro-sized dipole or magnetic trap ar- 
rays, we find for n = 70 a separation of the n = 2 
and the = 3 resonance peaks in the kilohertz regime, 
A2,3 = 27r X 146 kHz. Due to the strong scaling, even 
larger separations can be achieved for smaller lattice con- 
stants: at a = 5 /im, A2,3 = 27r x 9.4 MHz. These detun- 
ings need to be compared with the finite linewidth of the 
Rydberg state. For the n5i/2 Rydberg states of ^^Rb one 
finds Tn = 27r X 0.699 (n*)-^-^^ QHz where n* = n - (5 is 
the effective principal quantum number of the Rydberg 
state including the quantum defect 5 [21,. For our spe- 
cific example this yields Tjq = 27r x 3.0 kHz. Since the 
latter scales as and A^^^^^+i scales as n^^, the finite 
linewidth of the Rydberg state will eventually prevent the 
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resolution of the individual resonance peaks for smaller 
n. Nevertheless, there exists a broad parameter regime 
that does not put any restrictions on the experimental 
feasibility. 

The actual population of the desired many-body states 
that contain multiple Rydberg excitations relies on mul- 
tiphoton processes since - starting from the canonical 
ground state - Hamiltonian ([T]) can only couple many- 
body states that differ exactly by one Rydberg excita- 
tion. As a result, the excitation lasers need to be effec- 
tive long enough for a sufficient population of the desired 
final state. For the example of Fig. [2j we employed a 
maximal duration of T^ax = 301^"^. For n = 70 and a 
lattice constant of a = 5 /im this yields a time scale of 
3.6 /is which is more than one order of magnitude shorter 
than the corresponding Rydberg lifetime. Another pos- 
sible limitation of the proposed scheme stems from the 
uncertainty in the filling number A^o of the lattice. Any 
(site-dependent) variation in this number will affect the 
excitation dynamics by altering the collective Rabi fre- 
quency Q. However, for long enough excitation times T 
we do not expect a strong influence on our results since 
the excitation numbers and N^^ are predominantly 
determined by the detuning A of the laser rather than 
its coupling strength Q. We also note that Eq. ([6| is de- 
rived from Hamiltonian ([2| that neglects any influence 
of H\ cf. Eq. ([3|, on the spectrum. Including the latter 
slightly shifts [less than 2% in the case of Fig. [2] and splits 
apart the level crossings, cf. Fig. [2|c). As a result, the 
peaks of Figs. [2| a) and (b) develop a substructure which, 
however, can only be resolved for unrealistically long ex- 
citation times and a high resolution in A. If needed, 
the numerical determination of this substructure poses 
no difficulties. Because of this substructure, the widths 
of the resonances in Fig. |2] increase for increasing lattice 
sizes. 

While the present work focuses on the excitation of 
nSi/2 Rydberg states, a possible extension of the pro- 
posed scheme is the application to higher angular mo- 
mentum states for which the excitation cannot be mod- 
eled by means of a simple two-level system due to numer- 
ous possible molecular states with different symmetries. 
Instead of a single Cq coefficient one could possibly mea- 
sure an effective interaction strength that averages over 
all possible molecular symmetries and that is relevant 
for practical purposes [15 . On the other hand, different 
molecular interaction potential could also lead to sepa- 
rate peaks in the excitation spectrum which renders the 
extraction of a Cq coefficient ambiguous. 

In contrast, for the case of interacting nSi/2 Rydberg 
states one encounters a clean situation: The singlet and 
the triplet molecular states can be assumed degenerate 
for the interparticle distances we are considering [13 and 
any effects resulting from the hyperfine structure can be 
diminished by choosing the polarizations of the excita- 
tion lasers appropriately. For our example of ^^Rb, the 



excitation scheme 5Si/2{F = tuf = 2) ^ bPs/2{F = 
rriF =3) ^ nSi/2{F = tuf = 2) via circularly polar- 
ized lasers ensures that only these three hyperfine levels 
contribute to the excitation process. 

To conclude, we have shown that a coherently excited 
Rydberg gas in a one-dimensional lattice can be used for 
the determination of the interaction strength of Rydberg 
atoms. With the lattice spacing precisely known, the 
proposed scheme only relies on the measurement of the 
number of Rydberg atoms after an excitation cycle, cir- 
cumventing the mapping of the interaction energy as a 
function of distance. 

M.M. acknowledges financial support from the German 
Academic Exchange Service (DAAD). 



[1] T. F. Gallagher, Rydberg Atoms (Gambridge University 

Press, Gambridge, U.K., 1994). 
[2] D. Jaksch et a/., Phys. Rev. Lett. 85, 2208 (2000); M. D. 

Lukin et a/., ibid. 87, 037901 (2001). 
[3] D. Tong et al, Phys. Rev. Lett. 93, 063001 (2004); 

K. Singer et a/., ibid. 93, 163001 (2004); T. Gubel 

Liebisch et al, ibid. 95, 253002 (2005); T. Vogt et al, 

ibid. 99, 073002 (2007); G. S. E. van Ditzhuijzen et a/., 

ibid. 100, 243201 (2008). 
[4] R. Heidemann et a/., Phys. Rev. Lett. 99, 163601 (2007); 

M. Reetz-Lamour et a/., ibid. 100, 253001 (2008); T. A. 

Johnson et al, ibid. 100, 113003 (2008). 
[5] E. Urban et a/., Nat. Phys. 5, 110 (2009); A. Gaetan 

et al, ibid. 5, 115 (2009). 
[6] M. Saffman, T. G. Walker, and K. M0lmer, Rev. Mod. 

Phys. 82, 2313 (2010). 
[7] T. Wilk et al, Phys. Rev. Lett. 104, 010502 (2010); 

L. Isenhower et a/., ibid. 104, 010503 (2010). 
[8] B. Sun and F. Robicheaux, New J. Phys. 10, 045032 

(2008) 

[9] B. Olmos, R. Gonzalez-Ferez, and L Lesanovsky, Phys. 

Rev. A 79, 043419 (2009) 
[10] T. Pohl, E. Dernier, and M. D. Lukin, Phys. Rev. Lett. 

104, 043002 (2010). 
[11] B. Olmos, R. Gonzalez-Ferez, and L Lesanovsky, Phys. 

Rev. Lett. 103, 185302 (2009). 
[12] H. Weimer et a/., Nat. Phys. 6, 382 (2010). 
[13] K. Singer et a/., J. Phys. B 38, S295 (2005). 
[14] G. Boisseau, L Simbotin, and R. Gote, Phys. Rev. Lett. 

88, 133004 (2002); A. Reinhard et a/., Phys. Rev. A 75, 

032712 (2007); T. G. Walker and M. Saffman, ibid. 77, 

032723 (2008). 

[15] T. Amthor et a/., Phys. Rev. Lett. 98, 023004 (2007). 
[16] G. Ates et al, Phys. Rev. Lett. 98, 023002 (2007) 
[17] T. Amthor et al, Phys. Rev. Lett. 104, 013001 (2010) 
[18] N. Tezak, M. Mayle, and P. Schmelcher, 

arXiv :1012.3810y1 [physics . at om-ph] . 
[19] M. Mayle, L Lesanovsky, and P. Schmelcher, J. Phys. B 

43, 155003 (2010). 
[20] R. M. Potvliege and G. S. Adams, New J. Phys. 8, 163 

(2006) 

[21] A. L. de Oliveira et al, Phys. Rev. A 65, 031401(R) 
(2002). 



